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1. Da.scuss3.Qa 

Several recent experiments have shown features that indicate the presence 
of an H a ~ ion. Schultz 1 proposed a state of H 2 " in the formation of E and H“ 
by electron impact on hydrogen gas. Khvostenko and Duixel^hii 23 proposed a 

similar mechanism in electron impact measurements on H 2 , i.e., H 2 + e > (E 2 “) - 

3 H~. They also detected a negative ion of mass two by mass sp e ctroaet e r 

upon bombarding a mixture of water and antimony gas with electrons. In ref. 

(2b) i3 explained why this was felt to be H 2 “ with a lifetime long enough 
to be observed, rather than the other possibility, D“. In astrophysics 3 it 
lias been suggested that H 2 “ occurs in the solar corona and dissociates into K 
and H", 



The latest of these have been some electron scattering experiments by 
Simpson and Kuyatt, 4 * 5 where a highly monochromatic electron beam was scattered 
off hydrogen gas. A plot of current vs. energy was made, and resonances correa - 
ponding to a set of vibrational levels were observed* From these levels the 
spectroscopic constants ce e , m e x e , and the approximate Eo were calculated. Hie 

spectroscopic constants are very similar to those of H s , % or s T^j, ■ however the 
threshold energy was lower than both, thereby suggesting the possibility cf a 
resonant state of H 2 “ consisting roughly of an inner a electron and two outer 
electrons, with cynmetry %•* (Note that the outermost, planetary electron would 
not greatly affect the field in which the nuclei move.) 

The predominant common feature of all these experiments is that they 
are of a scattering nature, and thus the observations could be of a real bound 
excited state, or of a resonance in the continuum. 

Previous papers concerning theoretical calculations 6 on H 2 " show that 
there is no bound state with an energy lower than that of the H 2 X 1 £g state 
of the hydrogen molecule with an electron at infinity. In the computations 
of (6) many configurations were tried, with complete variation of non-linear 
and linear parameters, and in all cases the variation method itself produced 
the equivalent of a hydrogen molecule with an electron in the continuum. 

Thus it was decided to look for a bound state of a different nature, i.e., 
following the suggestions of experiment, a state with one inner a electron and 
two outer Tf electrons, such that essentially the wave function would resemble 
the trapping of an electron in a ir orbital by a 3 l^ i or hydrogen molecule. 

This appears physically resonable since an incoming electron would "see' 1 a o 
electron about /two nuclei, and a It r electron perpendicular to it, or therefore, 



* Z* and Z~ are essentially degenerate for this discussion. This can be seen 
from the small overlap of the (jj,q) parts of Itf and 3 ■%'„ 
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a net attractive charge which was greater thsn aero, uinra ;Ue exc-iteei. c r ,ho- 
gonal outer electron would art screen for the incoming one. Also,, in be tv 
the two electron (lo Iv) and the three election (lor $ir) systems, the tieotccns 
are highly uncorrelated. The importance cf this will he di 3 cussed late: 

Note the difference between this end electron capture by a hydro* 
atom to fom H~. In atoos an incoming electron can induce polarization art i.r 
the hydrogen molecule, most of the electron density laust regain between 'he t-’ 
nuclei to form the bond, ana therefore cannot be polarized nearly as nr. as 
two separate atoms by the incoming electron . This gives another physic: pic , .re 

of the reason a state like the 1 Eg state of hydrogen cannot capture an ?ctr 1 
but states like ^TTu and 3 T5[ 1 night be likely to. 

Use of a trial variation function of this type could lead to two ossi e 

results. First, but highly unlikely due to experimental evidence end p: riots 

theoretical calculations, there could be a true bound state which would e a 

solution to the Schroedinger equation, with a true variation principle v • l ar 

* 

energy less than that of the ground state of the hydrogen uolecule. Or ncor .y, 
the lo lif configuration could cause a narrow, but very deep potential w .1 Ter 
the incoming electron, such that there is a resonance, composed of a hi vj 
localized, relatively long-lived combination of continuum states (wave eket) 

The problem then is to construct such a highly localized wave par'- at 
from continuum functions having a lifetime long enough to he observable in 
other words, with a low probability of a radiationless transition from ;I e 
resonant state to the hydrogen molecule plus an elect ron of proper er.er ; in 
the continuum. 

Using the physical arguments above, a configuration which will hopefully 
have a deep potential well for the scattered electron and a smell overl i ■ wi h 
the Kg 1 F.g + e must be constructed. A preliml nary calculation was trade >ith r 
trial function with configuration loir hr . Even this rough calculation ’.cnf:. msd 
the suspicion that there is no bound discrete state with an energy less l;fcsn that 
of Ha 1 fg + e» The second possibility then becomes the one which, cur:- ; in- 
vestigated. Certain difficulties arise, for there is no longer s \arii ' on 
theorem and no rigorous minimum principle . How does one select the pas ete r-t 
of a wave packet and calculate the energy under these circumstances? 1 ar - - 
rents will be given to show the validity of an approximate variation pz Lucipi . 
They are based on the fact that the systes is to a large degree uneorrd ited , 
with weak interactions between electrons, and that the deep potential causes 
resonance packet, which, is t\ function of time, and this tine dependence ray, to 



&. first approxioation, oe «s-’pr-?3Sfcc :.n t;< ? jasie ,-s :f ha trial f ..a to.,« 

•'fel ther of tvies? arguments can be proven ..*1 joyously. but the submsquen Ib.s- 
easslon and res-olts show th.-ai both to be physically quite valid. 

The variation principle requires .'-.fes-p all bound states be ortV.c ; lal 
to all lower energy states sf the same s^im etry. In the j resent j-ase, he 
packet approximating a boon 1 state must b:-: orthogonal to nil s Eg -stile of 
H a ~Zg plus an electron at infinity. Aluheugb. rigorous to thogona Lii,; Jjr- 
possible since both states are effective orthogonal? l y (zex > ower 5 } ;• 

possible if the racket, is constructed so the two inner electrons xrt V 
wnict is orthogonal to Ik. 1 Lg. The variation is rental cht & such ;b 1 . ! 
inner electrons rust remin lag ltr u , i.e.» h a X 1T U and the: efore ?s • 0 

orthogonal to H 2 a £g. Orthogonality is also expected because tfc : cor.o 
ele ctron In the* lover state overlaps? negligibly with the 'ocaliz 1 Ity 
packet. TM 3 has been verified by direct c&lculat .on of he ave lap i 
of non-equal energy cases. If the packs'; vers to overlap with a *y c 
nteta it. would have to be those’ of equal energy. The fol owing . rgiv-c 
presented to account for thin case. It aut been notice d . a .retx jpei •' 
good agreement with experiment, the insensitivity of tba ariatk i cf 
trends parameters to the ether’s parameter?-, end the restating t -li »e 
molecular orbitals indicate -chat the ele :t:*an correlation is very seal 
tills is the case, l/r 1?i can be considers!, ns a perturbation and a zero •*r«le * 
argument that the overlap is zero can be /given. Higher order p« iturba on tebry 
will gi/e negligible overlap if there is ’’convergence" by zero order. f 
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e.nd the 



The packet is lag iir u x ,5in* « the cont.Lisr.um is la lr u ' 

clearly orthogonal by the orthogonality of He’"* orbitals. There! ore hj 
bat-ion theory, weak correlation implies tout a restricted variation '1 
to different H a * orbitals) allows a true adnimuen principle with which 
late the ’best" energy. To verify this iU is noted, in Eoj^eln, ur.d ?x 
Bailer's work ca helium 7 very reasonable results are obtained when Was 
principle is applied to such vn correlated packet states a *, (£p] 1 ' % Ig 
(2j-3p) l Fg, (2n^d) 1 D^ 1 , end {£p5> i) 3 B a „ the otter han . Hol^J or f:’r. i tii -.l 
the (2a) s ~ X S atate 7 ^ in H" hae a lover energy than that ov* the (/ s) 4 3 - stain 
cf 3~. In helium, the reverse > expected result is fonnu. It is fell 1. tne 
authors that the result on H™ is in error because in the .:•£») 2 ~S cess aorroU 5 . •• 
tion is not small and no approximate variation principle: Isolds. Clear: (Hr.) 1 
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being lower than (is) 2 is unphysieal and contrary to the more extensive calcu- 
lation of Pekerio® 2 where he shows only one bound state exists. The (2s) 2 1 S 



to the correlated nature of the system and the consequent failure of the 
variation method. 

Another Justification of the variation principle is that since this is 
a resonance, the energy calculated by variation is not a discrete eigenvalue, 
but the average of a spread of energies. This average energy must be a function 
of time, since the packet must eventually spread. Now there is no explicit 
time dependence in the variation function, but the variable parameters regulate 
the spread of the wave packet, and thus the parameters themselves may be con- 
sidered as time dependent and tending to cause the resonant wave packet to 
delocalize as time Increases. 

To achieve the longest lived resonant state in the restricted region 
of Hilbert space selected, therefore, by the uncertainty principle, AE av must 

a 

be minimized. 

This can be stated as E^ v EayCc^), where Oj_ are the variable parameters, 
and Q!j_ » a*(t). Then to minimize AS, 



But the variation method causes =* 0 for all therefore min- 

imizing AEav and maximizing At. This in essence selects the best (longest 
lifetime) wave packet possible in the region of Hilbert space considered. Note 
that the variation alone will not give a long lifetime state, but will only 
give the packet of longest lifetime among the class chosen. 

The use of the variation principle has now been physically Justified. 
Thus the question becomes: is the state which has been found of long enough 

lifetime to be observed? An appeal to the ultimate agreement with experiment 
could be made, with the claim that the state clearly exists long enough to 
vibrate. However, it is preferable to give a physical argument based on the 
weakness of the correlations and the resulting perturbation method to indicate 
a long lifetime. 

To demonstrate this, consider the standard auto-ionization picture of 
bound states in the continuum. 9 The argument will be outlined here for helium 
because it is easier; it will be clearly extendable to the case. If 



problem in helium is also not completely settled. 8a >^ This is felt to be; due 




(LI) 
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lA*i2 3 v corr. is small the zero order wave functions are products of two 
hydrogenic orbitals . The functions and levels fall into three classes with 
overlapping energies: 

(a) E = ~ ~a — --'2 both one -electron functions bound 

“a n 2 

1 ko*^ 

(b) E 1 + — one one-electron function bound and one continuum 

n x 2 

k 2 k» £ 

(c) E « both continuum 

When n x and n s are large in (a) and n x small in (b ) , a k can be found such 
that E con t, « Em,» i.e., there exists a bound state in the contiiiuum. A non- 
radiative transition from the bound to continuum state of equal energy is given 
by the probability of auto-ionization (see ref. 9, and ref. 7a). Therefore the 
transition probability, w, is proportional to 

j^onnd ( ~ | cont.)| 2 6 (Em. - ®eontJ 

The lifetime of the bound state is proportional to V 1 . 

Now if the correlation is small, an assumption upon which the equation 
for w was derived, the first order correction to the bound state function due to 
correlation is given by 

(Ho -W *1 * (^ - El) -0 

If Tjfi is expanded as usual on all other states of the system, and substi- 
tuted back into the perturbation equation, then multiplied by \|f contt and integrated, 
the result is,* 

^ bound j ^-| cont.) + Cm. ca (Em. - E cont J « 0 
l4i cont. 

Now when Em, * E coat. ‘ fclie matrix element is to first order zero and w is 
then of the font* x5(x), x => Em. - E cont. > aa<i therefore approximately zero. 

Higher order perturbation theory modifies this, but if the perturbation is really 
small, in a negligible way. 

In support of this argument it is pointed out that the types of resonant 
states that have been observed are generally highly un correlated and ones in 
which the "resonant trapped" electron sees a deep well due to the rest of the 
system. It is exactly this type of state for which successful variation 
calculations have been carried out. 



* In this expression, the continuum normalization factors have been cancelled 
from both sides. 
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II. Wave Function and Calculation Methods 



Using the physical arguments above, it was decided to use a one con- 
figuration three -electron variation function, having 2 £g symmetry, such that 
loss of the outer electron would not he able to yield a 1 S g state, but only 

The elliptic molecular orbitals and calculation methods of Harris 



a *17^, state 



and Taylor were used. 
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For one configuration the variation function becomes: 

Y <?,*)’ 6*i£Uj,(P)$<W €L-Co> 



(II. 1) 



where Uj^(P) is the Jd element of the matrix representation of the permutation 
P on the spin eigenfunction bases. The €l 



H 



are: 11 



(o' K - (&<*) 



The spatial part is the product, 

=> IT 0;&i) +TT 

where the are of the form. 



(II. 2) 



(II.3) 



r-ti 






and the 0^’ are, 

(- 0 0; ( f, - *s v\,vr>, -v ) 



(II. 5) 



T] and 0 being the standard two-center elliptical coordinates. The 5j, O-i 
are variable within the computer program used, but n,m,v are input parameters 
and therefore fixed throughout the problem. The configuration chosen was: 

°P,0 ♦ XO t - 0 -V 0' 0'0 ? ' (11*6) 



It is convenient to denote these in terms of the v±, since none of the n,m,v 
change and all » 0. Thus, 

Since Unj * 0 and £|vjj =2, there is no sign change between the products, 
and gerade symmetry is ensured. Also, 2V* = 0, so the three electron product 
is a £ state. 

Note that if the outer electron is removed (placed in the continuum), 

L&0, * . (H.8) 

is obtained, which is a IT state (£v ^ 0) with ungerade symmetry since the 0 1 ' 1 



7 



introduces a - sign between the products. Thus it has been spatially ensured 
that the 1 configuration variation function cannot vary itself to H 2 1 2g + e 
by varying 5-i to a very small number which has the effect of ejecting the 
outer electron into the continuum. 

Looting at the spin portion, it is also noted that ©^ 0 *. can "lose” an 
electron to become, 

0*,<UV> = (n.9) 

which is the spin eigenfunction for a two electron singlet. 

The one -configuration variation function is now constrained in such a 
manner that it cannot eject the outer electron and become anything but H s 1 1T U 
+ e, and the resonant wave packet is suitably localized by these constraints. 
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III . Calculations and Results 



In order to obtain a potential curve and the 
spectroscopic constants for comparison with those reported by experiment, 4 it 
>•73.3 planned to calculate five energy values a'o nuclear radii corresponding to 
scaled roots of the fifth order Chebyschev polynomial T 5 (x), where x = 

R - R e /. 34026, which gives separations of ± .2000 and ± o23'5 au. about R G , 
and fit these to Chebyschev polynomials To through T 4 . 1£ The spectroscopic 
constants cculd then be calculated by the method of Dunham. 13 

Preliminary calculations with a trial function of the type lo(Tw) 2 
yielded an equilibrium radius of 1.9196 au., and the final full variation of the 
5 and a parameters was done about this R e . Upon fitting to the Chebyschev 
polynomials and minimizing E(R) with respect to R, Rg was found to be 1.94-9 au. 
instead of 1.9196, so Dunham's method, based on a polynomial in (R - R^’f^, 
becomes excessively tedious. Therefore, the polynomial in x wa3 converted 
to a polynomial in R and the vibration-rotation constants calculated by the 
method used by Fraga and Ransil in ref. 14. 

To note the effect of greater interatomic distances three other R values 
were used. Table I gives the potential curve and parameters of the wave 
function. 

Table II is a comparison of the calculated spectroscopic constants with 
those reported in the scattering experiment. 

It is also of interest to note in Table III that as R increases, the 
potential curve approaches the curve of the X F U Gtate of H 2 ° 15 This can be 
explained by inspection of the parameters of the three electrons. At R * 10.0 au. 
there is essentially a Is hydrogen atom and a lp hydrogen atom widely separated, 
with on electron in between. In H 2 1 7T tl at this internuclear distance, the 
expected configuration is a pair of hydrogen atoms, one in the Is state and one 
in the lp state. 

This seems to indicate that H a ~ cannot be formed by collision of H and H“, 
but only by collision of Ha with an electron. Also, the proposed mechanism for 
formation of H and H“ by dissociation of Ha" may have another step, 

+ e «*-£> (H a ) 5 (Big + Bap + e) — £> H + H 

where the Hap nay recapture the electron of the Ha" or capture another of the 
electrons In the experimental system. 



9 



IV o Conclusions 

r 

The potential curve and electronic wave function for the resonant 2 Zg 
state of H 2 “ have been calculated by treating a one-configuration trial function 
as if it were a real variation function which could be minimised to a real 
discrete eigenvalue. An average energy very near to the experimentally ob- 
served energy was obtained, however, since this is a state in the continuum 
there ie no rigorous minimum principle, and the calculated average energy may 
be above or below the actual energy of the system. 

There is no doubt that if the configuration determining parameters n, 
m, and v were also allowed to vary, the final result would have been the 

one configuration approximation to the ground state energy of the hydrogen 
molecule. 

Excellent agreement with experimental spectroscopic constants has been 
obtained, which substantiates the statement that this is a resonance with a 
lifetime long enough to be observable, and the use of physical intuition to 
pick a restricted variation function which simulates an electron scattering 
off a deep, narrow potential well formed by an excited state of H 2 . 

The method of finding the longest -lived wave packet in the continuum 
by a restricted variation method coupled with good physical guesses may have 
considerable application in other problems. 

It is the intuitive feeling of the authors that all resonant states will 
be able to be treated by this approximate method, because if they have a long 
enough lifetime to be observed, they are probably weakly correlated, in the 
sense described above and therefore amenable to approximate variation treatment. 
This may not be a necessary condition, bat it is certainly a sufficient one. 

The lack of rigor prescribes great care in any application. 
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Table I a 

Wave Function and Potential Curve 



R (au) 


-S (au) 


8 


a 


R (au) 


-E (au) 


s 


a 


l«59o 


.74336 


1.115 

.445 

.100 


.765 

-.17 

0.0 


2.2432 


.74488 


1.420 

.635 

.105 


1.020 

--255 

0,0 


1.7196 


.74848 


1.175 

.485 

.105 


.810 

-.185 

0.0 


3.000 


.71513 


1.75 

.82 

.11 


1.34 

-.38 

0.0 


1.9196 


.75064 


1.270 

.550 

.105 


.890 

-.210 

0.0 


7.000 


.63377 


5.53 

1.78 

.11 


3.39 

-1.66 

0.0 


2.1196 


.74804 


1.365 

.600 

.105 


.970 

-.240 

0.0 


10.000 


.62929 


5.01 

2.51 

.11 


4.96 

-2.50 

0.0 










i5.oo b 


.61352 


6.24 

3.76 

.10 


6.21 

-3.76 

0.0 



a. Throughout this table, the electron parameters are listed in 
the order / 0 0 0\ so that the inner electron is first, and 

001 
\ 2 0 - 1 / 

the outer resonant electron is last. 

b. This point was not completely varied, but a definite continuation 
of the trend to Hi s + Hip + e“ is exhibited. 



Table II 

Comparison of Spectroscopic Constants 





calculated 


experimental 


E 


-o 751^8 au 


-.738 au 




lo9^9 au 




®e 


oJO ev 


.31 ev 




,042 ev 


a'.OI ev 


B e 


o0059 ev 




«e 


3 x 10“ 4 ev 




— £g- 


5->36 ev 





Table III 

Comparison with H 2 1 TT U Potential Curve 



R (au) 


-E H 2 C a 17 u a 


-E H 2 ‘ 2 Z g 


AS 


1.9 


.71422 


-75057 


.03635 


2.0 


.71418 


.74800 


.03382 


2.25 


.71006 


.74481 


.03475 


3.0 


.68445 


.71513 


.03068 


7.0 


.62478 


.63377 


.00899 


10.0 


.62448 


.62529 


.00481 



a. From ret, (12) 
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